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ABSTRACT 

The highly irradiated transiting exoplanet, HAT-P-7b, currently provides one of the 
best opportunities for studying planetary emission in the optical and infrared wave- 
lengths. We observe six near-consecutive secondary eclipses of HAT-P-7b at optical 
wavelengths with the EPOXI spacecraft. We place an upper limit on the relative 
eclipse depth of 0.055% (95% confidence). We also analyze Spitzer observations of 
the same target in the infrared, obtaining secondary eclipse depths of 0.098 ± 0.017%, 
0.159 ± 0.022%, 0.245 ± 0.031% and 0.225 ± 0.052% in the 3.6, 4.5, 5.8 and 8.0 micron 
IRAC bands respectively. We combine these measurements with the recently published 
Kepler secondary eclipse measurement, and generate atmospheric models for the day- 
side of the planet that are consistent with both the optical and infrared measurements. 
The data are best fit by models with a temperature inversion, as expected from the high 
incident flux. The models predict a low optical albedo of < 0.13, with subsolar abun- 
dances of Na, K, TiO and VO. We also find that the best fitting models predict that 
10% of the absorbed stellar fiux is redistributed to the night side of the planet, which 
is qualitatively consistent with the inefficient day-night redistribution apparent in the 
Kepler phase curve. Models without thermal inversions fit the data only at the 1.25-0" 
level, and also require an overabundance of methane, which is not expected in the very 
hot atmosphere of HAT-P-7b. We also analyze the eight transits of HAT-P-7b present 
in the E'POX/dataset and improve the constraints on the system parameters, finding a 
period of P = 2.2047308 ± 0.0000025 days, a stehar radius of = 1.824 ± O.O89i?0, a 
planetary radius of Rp = 1.342 it 0.068i?jup and an inclination of i = 85.7l2'2 deg. 
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Introduction 



The characterization of exoplanetary atmospheres is a quickly maturing field. Thus far, ob- 
servations have revealed an intriguing range of results that are in turn feeding back into a rapidly 
evolving theoretical paradigm. The most fruitful method of investigation to date is the detec- 
tion of secondary eclipses of hot Jupiters — highly irradiated gas giant planets in short period (2-3 
day) orbits. By measuring the relative depths of the secondary eclipse of an exoplanet in multiple 
bandpasses, one can reconstruct a low-resolution emission spectrum from the dayside of the planet. 

The planet/star flux ratio is highest at infrared wavelengths, due to the thermal emission 
from the planet. The Spitzer Space Telescope has been used to measure the infrared secondary 
eclipse depths for a growing list of extrasolar planets. These measurements reveal that the atmo- 
spheres span two distinct categori es: thos e with temperature in versions in the upper atmosph eres 
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It has been proposed that the distinguishing factor between the two 
classes is the incident flux on the planet, with higher incident flux resulting in gas-phase TiO /VO 
([Fortnev et al.l 120081 : iBurrows et al.l 120081) or non-equilibrium pr oducts of photochemistry, as seen 
in solar system bodies (jBurrows et al.l l2008l : IZahnle et al.ll2009l ). The additional absorber in the 
upper atmosphere generates the temperature inversion. Orbiting an F6 star, with a large stel- 
lar radius and high stellar temperature, in a short orbital period, HAT-P-7b is one of the most 
highly irradiated hot Jupiters, receiving 4.5 x 109 erg s ^ cm It is therefore predicted to have a 
temperature inversion. 

At optical wavelengths, the planet/star flux ratio is an order of magnitude less favorable than 
in the infrared, due to the much lower thermal emission from the planet at the shorter wavelengths, 
and also the low optical reflectance expected from cloud-free atmosphere models. Although ob- 
serving at optical wavelengths makes the detection of the secondary eclipse more difficult, it can 
provide very useful information. For instance, in the case of HD 209458b, the observed tempera- 
ture in version implies the p resence of an unk nown absorber hig h in the atmosphere, as mentioned 
above (jBurrows et al.l 120071 ). Subsequently, iRowe et al.l (j2008l ) placed very stringent upper lim- 
its on the optical depth of the secondary eclipse using the Microvariablity and Oscillations of 
Stars {MOST) satellite, and concluded that the absorber must have low reflectance. Optical sec- 



onda ry echpses have been detected for CoRoT-lb (jSnellen et al.ll2009l ) and CoRoT-2b (jAlonso et al 



20091 ) with relatively low signifi cance using the CoRoT spacecraft, and with high significance for 
HAT-P-7b (jBorucki et al.ll2009l ) using the iTep/er spacecraft. The only detection of an optical sec- 



ondary eclipse from the gro und at the time of submission is for OGLE-TR-56b in the z' band 



(ISing k L6Dez-Moralesll2009l ). 



The next step is to combine the optical and infrared secondary eclipse measurements to extend 
the planetary emission spectrum and refine statements about the atmospheric properties. Prior to 
the launch of Kepler, observations of HAT-P-7b were obtained with the EPOCh (Extrasolar Planet 
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Observation and Characterization) component of the NASA EPOXI Mission of Opportunity. One 
of the science goals of this mission was to measure or place meaningful upper limits upon the 
geometric albedo of a small set of transiting extrasolar planets. The goal of this paper is to 
combine the HAT-P-7b secondary eclipse measurements made with the EPOXI, Kepler and Spitzer 
spacecraft to generate a broadband spectrum covering 0.35 — 8/im, and to identify a set of model 
planetary atmospheres that reproduce the observations. The paper is organized as follows: In 
Section [21 we describe the EPOXI observations and data analysis; in Section [3l we describe our 
analysis of the Spitzer observations of HAT-P-7; in Section H] we combine the results of these two 
data sets with the published Kepler secondary eclipse measurement and describe the model fits to 
the final spectrum; and in Section [5] we discuss the conclusions of the paper and the extension of 
this analysis to additional EPOXI targets. 



2. EPOXI observations and analysis 

HAT-P-7b was the final transiting extrasolar planet target observed with EPOXI. We observed 
eight transits: the first in the 'pre-look', an initial two day observation window (UT 2008 August 
8 and 9) to refine the pointing, the latter seven observed consecutively in a fifteen day observation 
run lasting from UT 2008 August 16 to August 31. Seven secondary eclipses were observed during 
the run. We obtained 26,781 observations of HAT-P-7 in total. 

We observed HAT-P-7 with the 30-cm diameter High Resolution Instrument (HRI) on-board 
the OX/ spacecraft. We used the 1024x1024 pixel CCD with a clear filter, providing wavelength 
coverage from 350 to 1000 nm. The instrument is significantly defocused, with the PSF having a 
FWHM of 9.0 pixels (3.7 arcsec). Due to on-board memory limiting the number of full frame images 
that could be stored, we observed in two sub-array modes: 128x128 pixels typically, and 256x256 
pixels during times of select transits and eclipses to minimize pointing losses. Observations were 
taken nearly continuously, with gaps of approximately 7 hours every 2-3 days to download the data 
from the spacecraft. 



The data reduction pipeline is detailed in iBallard et al.l ((20091) and summarized here. We 



receive calibrated data products from the existing Deep Impact pipeline (jKlaasen et al.ll2005l ). to 
which we apply an additional bias correction. We use PSF-fitting to find the position of the target 
star in each image to within a hundredth of a pixel; we discard 4850 images at this stage due to the 
star falling out of the field of view or too close to the edge of the CCD for accurate photometry. We 
discard an additional 40 images due to significant (15-cj) deviations from the model PSF, caused 
by readout smear or radiation events. We then apply several corrections based on our knowledge 
of the CCD architecture, including scaling down the fiux values in the two rows and two columns 
that form the internal boundaries at the center of the CCD, which has readout electronics in each 
of the four corners. We also scale the counts in images taken in 256x256 mode to those taken 
in 128x128 mode, to account for a systematic fiux offset between the two modes. Finally, we 
divide the images by a fiat field generated from a green LED stimulator lamp mounted beside the 
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CCD, which corrects for changes in the position-dependent sensitivity that have occurred since 
the detector was characterized prior to launch. We obtained these calibration frames ('stims') in 
six short blocks, each of approximately 20 minutes duration, spaced every few days during the 
observations of HAT-P-7. Each block contained 250 images, and in total we obtained 500 stims in 
128x128 mode and 1000 in 256x256 mode. We correct each science frame in a given mode with 
the stims obtained in that mode. 

We perform aperture photometry on the resulting images, using a 10-pixel radius aperture. At 
this stage there is still considerable correlated noise in the data, shown as the lower light curve in 
Figure[Tl due to the time- variable jitter in the spacecraft pointing. For the observations of HAT-P-7 
this jitter was on the order of 10 arcsec (40 pixels) hr~^. We find that the monochromatic green 
stims are a poor representation of the interpixel sensitivity, which seems to have a color-dependent 
component, and as a result the pointing jitter introduces correlated noise that we are unable to 
remove with the calibration frames. We subsequently use the out-of-transit and out-of-eclipse data 
to map out the brightness variations with CCD position. For a small random subset (6% of the 
21684 points in the light curve) we find an iteratively sigma-clipped mean of the closest spatially 
neighbouring points, and use these values to generate a surface. We then fit a 2-dimensional spatial 
spline to the surface, and correct each point in the light curve by interpolating onto this surface. 
We discard an additional 207 frames affected by readout smear or radiation events at this stage, 
resulting in a final set of 21684 images. The full corrected light curve is shown as the upper light 
curve in Figure[Tl with an out-of-transit scatter that is 91% above that expected from Poisson noise. 
The excess is due to the remaining correlated noise that we have been unabl e to calibrate, and is 



higher than the 56% value we attained for the ii'POX/ observations of GJ 436 ([Ballard et al.ll2009l ). 
This is due to the fact that the HAT-P-7 light curve has approximately two-thirds the number of 
observations of GJ 436, and that these observations are distributed over a larger fraction of the 
CCD. Both of these effects decrease the robustness of the spatial spline and result in higher residual 
correlated noise. 

There are two 'fiare'-type events visible in the light curve, both of which we believe to be 
instrumental artifacts. The first occurs during the pre-look, at which time one of the four CCD 
quadrants was affected by a particular systematic flux correlation that was observed for other 
EPOXI targets. For the remainder of the HAT-P-7 observations this systematic effect was not 
present, and therefore the data obtained during the pre-look in this particular quadrant are not 
well calibrated by the spline. The second, which can be seen at approximately 18.6 days in Figuredl 
occurs at a time when the star was illuminating a part of the CCD to which it does not return in 
the course of the HAT-P-7 observations. The spline procedure we describe above therefore fails, 
since it relies on multiple visits to the same parts of the detector. 

In order to determine the limb-darkening coefficients across the -BP OX/ bandpass, we use the 
four stellar model atmospheres (Kurucz 1994, 2005) that bracket th e HAT-P-7 valu es of Tgfj = 



6350K and log g = 4.07, with [M/H] = 0.3 and vturb = 2.0 km/s (|Pal et al.l 120081 1. We found 



previously that for the wide EPOXI bandpass the choice of limb-darkening treatment did not 
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Fig. 1. — Upper panel: The full HAT-P-7 EPOXI W^t curve. The first transit was observed during 
the two day 'pre-look' that established the pointing; a further seven were observed consecutively 
over the following 15 day observation run. The lower light curve has not had the surface spline 
correction applied; the upper light curve includes the correction and is the final light curve used in 
the analysis. The two brightening events are instrumental in origin (see text). The vertical dashed 
lines indicate the expected secondary eclipses. Lower panel: The scatter in the out-of-transit data 
(shown as diamonds) does not decrease with increasing bin size as expected for Gaussian noise 
(1/VW, where N is the number of points in the bin, shown as the solid line normalized to the 
unbinned value of the scatter), indicating that there is a component of correlated noise remaining 
in the light curve. 



systematically affect t he resu l ts (IBallard et al.ll2009l ). We fit the four coefficients of the non-linear 
limb-darkening law of Claret (2000) to 17 positions across the stellar limb. We repeat this fit in 
0.2 nm intervals across the 350-1000 nm bandpass, weighted for the total sensitivity (including 
filter, optics and CCD response) and photon count at each wavelength. We calculate the final 
set of coefficients as the average of the weighted sum across the bandpass, bi- linearly interpolated 
between the four models, for which we find ci = 0.43, C2 = 0.14, C3 = 0.46, and C4 = —0.25. 



2.1. Transits 



We first consider the transits observed by EPOXI. We fit each transit with a linear, time- 
dependent slope using four hours of data before and four hours of data afte r the transit, to r emove 
any residual instrumental linear trend. We use the analytic algorithms of iMandel &: Agoll (j2002l ) 
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to fit the transit parameters, using as the goodness-of-fit estimator. Using the Levenberg- 
Marquardt algorithm, we fit for three geometric parameters: Rp/Ri,, where Rp and R^, are the 
planetary and stellar radii respectively; R-^/a, where a is the semi-major axis; and cos i, where 
i is the inclination of the orbit. We also fit for eight times of transit center for a total of eleven 
free parameters. We initially fix the period to the value from the discovery paper (jPal et al 
20081 ). 2.2047299 days. We find Rp/R^ = 0.0756 ± 0.0013, R^/a = 0.2249 ± 0.0116 and cos i = 



0.0754^QQgQg. In order to account for the effects of the remaining syster natics in the e rror b udget, 
we use the "rosary bead" residual permutation method described by I Winn et al.l (j2008l ). with 
2000 permutations. The quoted 1-a errors encompass 68% of the permutations. The light curve 
is shown in Figure [2l phased and binned in 5 minu te interval s , with the best fit model transit 
overplotted. Using the stellar mass value derived bv |Pal et al.l (120081 ). = 1.47 ± O.OSM©, we 
find R^ = 1.824 ± 0.089/?©, Rp = 1.342 ± 0.068i?jup and i = 85.7t|^ deg, which are consistent 
with and an improvement over the values in the discovery paper. However, they are m oderately 
inconsistent with the parameters presented in the recently posted IWinn et al.l (|2009l ) paper; a 
rigorous comparison of the results is not undertaken here. The new EPOXI values and the best fit 
times of transit center are given in Table [TJ For the transit times, shown in Figure [3l we see no 
evidence for transit timing variations on the order of 1-2 minutes, and we use a weighted linear 
fit to derive a new ephemeris combining the discovery ephemeris and the eight times listed here: 
rc(BJD) = 2454700.81308 ± 0.00100, P = 2.2047308 ± 0.0000025 days. 



2.2. Secondary eclipses 

We monitored eight secondary eclipses of HAT-P-7b with EPOXI. Two of the eclipses are 
interrupted by the systematic effects described earlier, and for the final analysis, we remove these 
two eclipses and use the remaining six, shown in Figure |H For this analysis, we perform an 
additional calibration step. When producing the full HAT-P-7 light curve, we use only 6% of the 
data points to create the 2-dimensional spatial spline, in order to minimize suppression of any 
additional astrophysical signals in the data, such as transits of additional planets in the system. 
However, when trying to extract the shallow signal of the eclipses of HAT-P-7b, we can make use 
of all of the data points that are not obtained during those events to increase the robustness of 
the calibration. For each event, we consider the points that occur within ±1.5 times the transit 
duration of HAT-P-7b from the predicted center of eclipse, resulting in a window of approximately 
12 hours duration. For each point in this window, we find a resistant mean of the fluxes of the 
nearest 30 data points on the CCD that do not occur during a transit or eclipse of HAT-P-7b, and 
normalize the flux of the point by this value. This reduced the scatter of the points in each window 
by 1-2% on average. Since we find that we are still ultimately limited by systematic errors in the 
resulting light curves, we use the scatter in the six measured depths as a proxy for our precision. 

For each eclipse, we bin the data in 20 minute intervals, rejecting 3-a outliers within each 
bin. As with the transit analysis, we fit the data on either side of the eclipse with a linear slope 
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Fig. 2. — The EPOXI HAT-F-7 light curve, phase-folded and binned in five minute intervals. The 
solid line shows the best fit transit model, and the lower panel shows the residuals after this model 
is subtracted from the data. 



with time to remove longer timescale trends. We generate a model using the geometric parameters 
derived from the transit fitting and with no limb-darkening, and scale the model depth to that 
which minimizes the of the fit to the binned data. We do not allow the time of secondary eclipse 
to vary, fixing it at the expected phase value of 0.5, since there is no evidence of eccentricity for 
HAT-P-7b (see Section[3]). We use the rosary bead method to determine the 1-a error bars, shifting 
the residuals and re-binning each time, with approximately 820 permutations per eclipse. The six 
eclipse depths are given in Table [21 where a negative eclipse depth corresponds to an increase in 
flux at the predicted time of secondary eclipse. The best fit eclipse models are overplotted in 
Figure HI The error bars in the figure are calculated from the scatter in each bin normalized by 
the square-root of the number of points in that bin. For the final value we take the mean of the 
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Fig. 3. — Observed transit ti nies of HAT-P-7b minus the best fit linear ephemeris. Upper panel: 
The sohd data point is from iPal et al.1 (|2008l ). the hollow points are the EPOXI transit times. 
Lower panel: An expanded view of the EPOXI transit times. 



six measurements; the systematics preclude us from scaling down the individual error bars when 
combining the measurements, resulting in a best fit depth of 0.0070 it 0.0240%, where the error is 
set such that 2-a (0.048%) encompasses the six points. 



3. Spitzer observations and analysis 



Secondary eclipses of HAT-P-7b were observed with the Spitzer Space Telescope (IWerner et al 



20041 ) under a Target of Opportunity program (program 40685) led by Joseph Harringto n. Observa- 



tions of two consecutive eclipses were obtained with the InfraRed Array Camera (IRAC: lFazio et al 
20041 ): the first was observed on UT 2008 October 28 in the 3.6 and 5.8 micron channels for 9.1 
hours, the second on UT 2008 October 30 in the 4.5 and 8.0 micron channels for 8.1 hours. The 
latter were preceded by a 30 minute pre flash set of observations, d e scribed further i n Sec t ion 13.21 
We treat the da t a in the sarn e fashion as ICharbonneau et al.l (j2008l ) ; iKnutson et al.l (j2008l . l2009al ) ; 
Machalek et al.l (j2008l . l2009l ). We summarize the procedure here, with particular care taken to 
describe any ways in which our analysis differs from earlier work. 

For each channel, we perform aperture photometry on the BCD (Basic Calibrated Data) prod- 
ucts, using an aperture radius of 2.0 pixels for the 3.6 and 5.8 micron channels, 2.5 pixels for the 
4.5 micron channel, and 3.0 pixels for the 8.0 micron channel. We determine the centers of the 
apertures by calculating the weighted sum of the flux in a circle of radius 3 pixels positioned at 
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Table 1. HAT-P-7 system parameters 



Parameter Value 



Adopted values^ 

M, (Mq) 1.47 ±0.08 

Mp (Mjup) 1.776 ± 0.077 



Derived values 
P (days) 
(BJD) 

R*iRe) 
RpiR jup) 
i (deg) 

Transit times (BJD) 



2.2047308 ± 0.0000025 
2, 454, 700.81308 ± 0.00100 
1.824 ±0.089 
1.342 ±0.068 
85.7+1:1 
2,454,687.5815 ±0.0012 
2, 454, 696.40413 ± 0.00075 
2,454,698.6064 ±0.0012 
2,454,700.8139 ±0.0011 
2, 454, 703.01513 ± 0.00099 
2, 454, 705.22219 ± 0.00086 
2, 454, 707.43077 ± 0.00074 
2, 454, 709.63062 ± 0.00092 



''Masses are from lPal et~al] ll2008h . 



Table 2. HAT-P-7b ii'P OX/ secondary eclipse measurements. Negative values imply increases in 

flux at the predicted time of eclipse. 



Predicted time of eclipse (BJD) 


Eclipse depth 


2,454,695.30127 


0.018 ± 0.014% 


2,454,697.50586 


-0.041 ±0.011% 


2,454,699.71045 


0.000 ± 0.015% 


2,454,701.91553 


0.016 ± 0.018% 


2,454,704.12026 


0.023 ± 0.010% 


2,454,708.52979 


0.025 ± 0.006% 


Average depth 


0.0070 ± 0.0240% 
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Fig. 4. — The relative flux versus orbital phase for the six HAT-P-7b eclipses observed with EPOXI, 
binned in 20 minute intervals, and offset by a constant for clarity. The solid line shows the best fit 
eclipse model in each case. 



the approximate location of the star. We choose the aperture size to minimize the scatter in the 
out-of-eclipse data. From this point, we treat the two shorter wavelength channels in a different 
manner to the two longer wavelengths and the analyses are described separately. 



3.1. 3.6 and 4.5 micron channels 



For each image, the position of the aperture was compared to the median of the preceding 
10 and following 10 images. Where the position deviated by more than 3 times the rms of these 
values from the median, the estimate of the centroid was corrupted and the image discarded. For 
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the 3.6 micron channel, 26 (0.53%) of 4906 images were discarded; for the 4.5 micron channel, 11 
(0.25%) of 4370 images were discarded. For both channels, the sky background was calculated in 
an annulus centered on the aperture with inner and outer radii of 10 and 15 pixels respectively and 
subtracted. 

Observatio ns in these shorter wavelength channels are characterized by a position- dependent 
flux sensitivity ([Reach et al.ll2005l : ICharbonneau et al.ll2005l : iKnutson et al.ll2008l . l2009al ): we fit for 
this sensitivity and the eclipse light cur v e sirn ultaneously. We generate the light curve with the 
analytic equations from iMandel &: Agoll (|2002l ) with no limb darkening, using the revised HAT- 
P-7b parameters from the EPOXI analysis. We model the intra-pixel flux sensitivity for both 
channels with linear and quadratic coefficients in x and y position, as shown in Equation [U where 
the measured flux /' is a function of the incident flux / and position; adding higher order position 
terms or time-dependent terms did not significantly improve the final fit. 



/' = f{hi + h2{x -x) + bs{x - xf + 64(y -y) + h{y - yf) (i) 

For each channel, there are seven free parameters in the fit: a constant, a linear and a quadratic 
coefficient in x position, a linear and quadratic coefficient in y position, an eclipse depth, and a 
time of ceii t er of eclipse. For th e fitting, we use the Monte Carlo Markov Chain (MCMC) method 



(|Fordll2005l : iHolman et al.ll2006l ): in each case we run four chains of 500,000 steps, varying one of 
the five parameters at each step, and discarding the first 20% of each chain. The best fit value 
of each parameter is taken as the median of the distribution, and the error as the range centered 
on the median that encompasses 68% of the total points. We find depths of 0.098 ± 0.017% and 
0.159 lb 0.022% in the 3.6 and 4.5 micron bands respectively. The results are shown in Figures O and 
El the former shows the uncorrected data in 6 minute bins, with the model of intra-pixel sensitivity 
overplotted, the latter the binned data with the sensitivity divided out, and the best fit eclipse 
model overplotted. In order to assess the effect of any remaining correlated noise in the data we 
also use the rosary bead method described in Section [2.11 In each case the number of permutations 
of the residuals is set to the number of points in the light curve (4880 points and 4359 points in the 
3.6 and 4.5 micron bands respectively) and the eclipse fit performed with the seven free parameters 
described above. We measure nearly identical eclipse depths and errors for the rosary bead method 
as for the MCMC method, finding 0.099 ±0.014% in the 3.6 micron band and 0.162 ±0.016% in the 
4.5 micron band. For the remainder of the paper we adopt the MCMC eclipse depths and errors. 



3.2. 5.8 and 8.0 micron channels 

The 5.8 and 8.0 micron channels show a long timescale time-dependent (as compared to 
position-dependent) flux sensitivity. The correlation is strongest in the 8.0 micron channel, where 
an increase in the flux during the course of the observations, referred to as the 'detector ramp' is 
evident. One interpretation of this is an increase in the effective gain of the pixels with time in 
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the 8.0 micron detecto r. One successful method for mitigating this ramp is the preflash technique 
(jKnutson et al.ll2009bl ). where a bright diffuse region is observed immediately prior to the start of 
the target observations in order to pre-emptively increase the effective gain of the pixels. For the 
8.0 micron observations of HAT-P-7, a 30 minute set of preflash observations, 3648 images in total, 
were obtained of the HII region GAL 075.83+00.40 (a = 20 21 39.5, 6 = +37 31 04). The median 
flux in the pixels in the aperture used for the HAT-P-7 observations ranges from 3500-4500 MJy 
sr~l. After the pre- flash a time-dependent correlation is evident in the data, albeit significantly 
attenuated when com pared with sim ilar time series analy s is per forme d in this channel (for example, 
Deming et al.l (j2007l )). As noted bv ICharbonneau et al.l (j2008l ) and iKnutson et al.l (j2008l ) we find 
that removing the most significantly affected data, in this case the first 18 minutes (79 of 2185 
images, 3.6% of the total), removes the need for higher order terms in the fit to the correlation 
described below. We also perform an iterative 3-fT fiux cut over the running median to discard 
outliers, possibly caused by radiation events, which removes 47 (1.9%) of 2453 images in the 5.8 
micron channel, and 51 (2.4%) of the remaining 2106 images in the 8.0 micron channels. 



We chose the time-dependent function that minimized the out-of-eclipse scatter, which was 
a linear function of ln(dt) in both channels, where dt is the time elapsed from the center of the 
first o bserved transit, shown in Equation [H Previous analyses have used quadratic fu nctions of 



ln{dt) (jCharbonneau et al 



2008 



Knutson et al 



20081 . 



2009a 



Machalek et al.ll2008l . l2009l ^. however 



by comparison the ramp is attenuated in both channels in these data, and adding higher order 
functions of time or position-dependent functions did not improve the fit. 



/' = /(ci + C2(ln(dt + 0.02))) (2) 



The additional factor of 0.02 prevents the singularity at dt = 0, and similarly to lKnutson et al 



(|2009al ) we find no significant change in the measured depth of the eclipse on varying this value. 
The MCMC fit was performed as described previously, with three free parameters in each case: a 
constant, a coefficient in time, and an eclipse depth. We found that when the time of eclipse was 
included as a free parameter there was a correlation between the time of eclipse and the coefficient 
in time. Since the measured times of eclipses in the shorter wavelengths do not deviate significantly 
from the predicted times for e = 0, we do not allow the time of eclipse to vary in the final fit for the 
5.8 and 8.0 micron channels. We instead fix them to the values predicted by the E'POX/ephemeris 
for zero eccentricity. We calculate the values and errors in the same fashion as for the 3.6 and 4.5 
micron channels, finding 0.245 ± 0.031% and 0.262 ± 0.027% in the 5.8 and 8.0 micron channels 
respectively. In order to confirm the validity of fixing the time of eclipse, we repeated the analysis, 
this time drawing the time of eclipse from a Gaussian distribution defined by the times of eclipses 
and errors measured in the simultaneously observed channels, using the 3.6 micron values for the 5.8 
micron channel and the 4.5 micron values for the 8.0 micron channel. We find near-identical results 
using this analysis. The results are shown in Figures [S] and O Again, we performed a rosary bead 
analysis to measure the effect of remaining correlated noise on the eclipse parameters. We find that 
a depth of 0.256 ± 0.029% in the 5.8 micron band, consistent with the depth measuring using the 
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MCMC method. In the 8.0 micron band the depth is 0.225 it 0.052%, approximately l-cr shallower 
and with a larger error than that produced by the MCMC method, indicating that correlated noise 
represents the dominant contribution to the error budget for this channel. Therefore we adopt the 
rosary bead method eclipse depth and error for the 8.0 micron band for the following analysis. 



1.00 



V) 

c 
o 
o 



> 



0.98 



I 0.97 



0.96 



~i — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 
3.6 micron 




4.5 micron 



5.8 micron 




5.0 micron 



0.95 Lj I I L 



J I I I I I I I L 




J I I L 



0.40 



0.45 



0.50 
Orbital Phase 



0.55 



0.60 



Fig. 5. — Secondary eclipses of HAT-P-7b observed in the four Spitzer IRAC bands. The data are 
binned in 6 minute intervals and offset in flux for clarity. The solid line is the best fit model of the 
position- and time-dependent instrument effects in each case, described in Section [3l 



4. Discussion 



Figure [7| shows the four Spitzer secondary eclipse depth measurements, the 95% confidence 
EPOXI upper limit at 0.65 micron, and th e recently published Kepler eclipse depth of 0.0130 ± 
0.0011% at 0.63 micron ([Borucki et al.ll2009l ). This latter measurement was derived from the initial 
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Fig. 6. — The same data as in FigureEl with the position- and time-dependent correlations removed. 
Overplotted is the best fit echpse model from the MCMC analysis in each case. The best fit depths 
and times of eclipse are given in Table [3l 



10 days of Kepler calibration data and is expected to be refined in the future. We use these eclipse 
depths to calculate the brightness temperature in each bandpass. Assuming the planet emits as a 
blackbody, this is a measure of the dayside temperature of the planet required to generate sufficient 
flux in the bandpass to reproduce the measured planet-to-star ffux ratio. In each bandpass, we 
integrate the ratio of a black body spectrum to a model stellar atmosphere o f HAT -P-7, with 
Teff = 6500K, log g = 4.0, [M/H] = 0.3 and wturb = 2.0 km/s (jKuruczl flQoi . bood l. weighted 
for the instrument throughput, and solve for the temperature required to match the measured 
eclipse depth. As expected, the highly irradiated HAT-P-7b has the hottest dayside brightness 
temperatures measured for any exoplanet, with temperatures in the four bands ranging from 2250K 
to 3190K. These values are given in Table [3l We also calculate this value for the Kepler secondary 
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eclipse depth measurement (jBorucki et al.ll2009l ) and find Tb = 3175 ± 40K 



We fit tfie Kepler and Spitzer data using tlie model described in iMadhusudhan &: Seager 



(|2009al ) but extended to vis ible wavelength s (Madhusudhan & Seager, in prep.), using a photo- 
sphere radius of 1.363i?jup (jPal et alJ 120081 ) . The model assumes global energy balance between 
the incident and emergent radiation. The HAT-P-7b atmosphere is too hot for condensates; the 
visible-wavelength opacities include TiO and VO, Na and K, and Rayleigh scattering. We include 
additional molecular opacities of H2O, CO, CH4, and H2-H2 collision-induced opacities (CO2 and 
NH3 are included in the model but are not required by th e current data). Our H2O, CH4, CO and 
NH3 molecular line data are from (jFreedman et al.ll2008l. and referenc e s ther ein). Our CO2 data 
are from R. Freedman (personal com munication) and Rothman et al. J 20051), an d we obtain the 
H2-H2 collision- induced opacities from Borvsow et al. ( 1997I ). and Borvsow ( 2002 ). 



Our model has fifteen free parameters, and we have at hand only five observational constraints. 
Consequently, the data allow a large number of solutions at any level of fit. Therefore, we nominally 
define 'best fit' models as those that fit the data to within the 1-cr errors shown in Figure [71 An 
exhaustive exploration of the fifteen dimensional model parame ter space is beyond the scope o f 
the present study. Nevertheless, based on strategies outlined in IMadhusudhan &: Seageii (j2009al ). 
we have been able to empirically assess the region in the parameter space that best explains the 
data. Figure [7] shows one representative model that provides a best fit to the observations. The 
red circles with error bars are the data and the green circles are the bandpass integrated model 
points. The red dotted lines show the Kepler and Spitzer bandpasses. The black curve is the 
model fit. The blue dashed curves show two black-body spectra corresponding to 2029 K and 
2974 K. The orange dashed curve shows a a black-body spectrum at 2600 K, for reference. The 
inset shows a zoom-in of the Kepler point and the model in the 0.4-1.0 fiui range. The cyan and 
purple curves in the inset show the thermal emission and scattered light, respectively. Figure [8] 
shows the pressure-temperature (P-T) profile of the model atmosphere. 

We find that all of our best fit models require a thermal inversion to fit the data. A thermal 
inversion is strongly suggested by the high flux ratio in the 4.5 micron channel of Spitzer compared 
to the 3.6 micron channel, in addition to the high flux ratios in the 5.8 and 8 micron channels. The 
emission features in the 4.5 micron and 5.8 micron channels are due to CO and H2O, respectively. 
The features in the 3.6 micron and 8 micron channels are predominantly due to CH4. The molecular 



Table 3. HAT-P-7b Spitzer IRAC secondary eclipse measurements 



Wavelength {fJ.m) 


Depth (MCMC) 


Depth (Rosary) 


Time (BJD) 


O-C (minutes) 


Tfc (K) 


3.6 


0.098 ± 0.017% 


0.099 ± 0.014% 


2,454, 768.05200 ± 0.0035 


7.20 ±5.04 


2250 ± 180 


4.5 


0.159 it 0.022% 


0.162 ±0.016% 


2,454, 770.26413 ± 0.0039 


-3.17 ±5.62 


2600 ± 210 


5.8 


0.245 ± 0.031% 


0.256 ±0.029% 






3190 ± 280 


8.0 


0.262 ±0.027% 


0.225 ±0.052% 






3100 ± 240 
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mixing ratios (with respect to H2) for the model shown in Figure[7l are H2O = 1.0~^, CO = 1.0~^, 
CH4 = 1.0^^, and no C02- If we relax the model assumption of thermochemical equilibrium, and 
allow a 'best fit' to the data to be defined by 1.25-ct errors, it is possible to generate suitable models 
that do not have a thermal inversion. However, these models require an overabundance of CH4 in 
the atmosphere. Given the very hot atmosphere of HAT-P-7b, thermochemical equilibrium favors 
a high abundance of CO and low CH4. The refore, a non-inversion rnodel fi t to data would imply 
signs of extreme non-equilibrium chemistry (jMadhusudhan &: Seagerll2009bl ). 



For a planet as hot as HAT-P-7, both scattered starlight and planet thermal emission could 
contribute to the planet flux in the Kepler bandpass (jLopez- Morales &: Seageri 120071 ). The relative 
amount of scattered starlight and thermal emission depends on the absorbers and scatterers at 
visible wavelengths: molecular opacities due to TiO and VO, atomic opacities due to Na and K, and 
Rayleigh scattering mostly from H2. HAT-P-7b is too hot for condensates, and without a reflective 
condensate layer, visible-wavelength photons will be absorbed before they can be scattered, except 
at very blue wavelengths. Even without a rigorous quantitative model, therefore, we can infer that 
the albedo of HAT-P-7b is likely to be low. 

Our best fit models show a low albedo and include the prominent sources of opacity in the 
optical. The model shown gives a geometric albedo of about 0.13 in the Kepler bandpass, and has 
scattered light dominating at wavelengths blueward of approximately 0.6 microns and the planet's 
thermal emission dominating redward of approximately 0.6 microns. The abundances of Na and K 
are free parameters in the model and are found to be 0.01 solar; similarly the TiO/VO abundances 
are found to be 10"'^ of their solar abundances. We note that even though TiO and VO are expected 
to contribute to t he visible and IR o pacity, whether they are responsible for the thermal inversion 
is not known (see lSpiegel et al.ll2009l ). 



We emphasize that the Kepler data point for HAT-P-7b was instrumental in constraining the 
P-T profile, beyond the constraints placed by the Spitzer data alone. This is particularly true 
because of the high planet-star flux ratio of HAT-P-7b; in general, a visible-wavelength observation 
with a lower contrast may not be as useful for constraining P-T proflle. 

The values of atomic and molecular abundances present in our models are valid in the frame- 
work of a ID averaged day-side atmosphere. The true values of atomic and molecular abundances, 
and hence the resulting albedo, in the atmosphere of HAT-P-7b depend on the 3D atmospheric 
structure, and in 3D are degenerate at the resolution of the Kepler data point. 



The i^Tep/er phase curve (jBorucki et al.ll2009l ) shows a very inefficient day- night redistribution 
of absorbed stellar energy. Although energy redistribution is not a focus of this paper, we emphasize 
that a relatively inefficient day-night redistribution is consistent with, and even required by, our 
best fit models. We find that approximately 10% of the incident stellar flux is redistributed to the 
night side, under the assumption that the visible and infrared thermal phase curves will be similar. 
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Fig. 7. — The EPOXI, Spitzer and Kepler secondary eclipse measurements of HAT-P-7b. The sohd 
Une is a representative best fit model, with a temperature inversion. The points with error bars 
are the observed measurements; the points without are the bandpass-integrated model values. The 
dashed lines show black-body spectra corresponding to 2029 K, 2600 K and 2974 K, for reference. 
The dotted lines are the instrument response curves. The inset is an expansion of the optical region 
of the spectrum, showing the Kepler measurement and response curve. The cyan and purple lines 
are the thermal and scattered components of the model respectively. 



5. Conclusions and future work 



We have assembled a broadband emission spectrum of the highly irradiated dayside of HAT- 
P-7b, extending from 0.35-8.0 micron. We find that this spectrum is best reproduced using an 
atmospheric model with a temperature inversion, supporting the hypothesis that the presence of 
an inversion layer is related to the absorbed stellar radiation. A census of exoplanet atmospheres 
at 3.6 and 4.5 micron that is being undertaken as part of the Warm Spitzer program will shed more 
light on the parameter space between where temperature inversions are present in the atmospheres 
of the most highly irradiated planets and absent in the less irradiated planets. This will include 
full phase curves of HAT-P-7 at 3.6 and 4.5 microns, which will further constrain the redistribution 
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Fig. 8. — The temperature-pressure profile of HAT-P-7b corresponding to the atmosphere model 
shown in Figure [71 showing a temperature inversion. 



and energy budget of the planetary atmosphere. In addition, as the Kepler observations of HAT- 
P-7 accumulate, the quality of the resulting light curve will be exquisite. This will not only 
provide higher constraints on the secondary eclipse and phase curve measurements, but will allow 
investigation into the variability of the planetary emission properties. 

At optical wavelengths, we find that the best fit models predict a low geometric albedo of 
< 0.13, which is far less reflective than gas giants in our solar system (e.g. 0.52 for Jupiter). The 
extremely hot atmosphere precludes the condensation of a reflective laye r (i.e. clouds), an d our 



predicted albedo is consistent with those measured for other hot Jupiters (jRowe et al.ll2008l ) 



The EPOXI HAT-P-7 secondary eclipse constraint is one of the few opportunities to confirm 
the Kepler measurement, given that the relative depth of the signal at visible wavelengths is on 
the order of 10^^. We observed a second transiting exoplanet target in the Kepler field of view 
with EPOXI, TrES-2, obtaining similar phase coverage to HAT-P-7, and will present those results 
in a forthcoming paper. Spitzer observations of TrES-2 have been obtained in the four IRAC 
bandpasses, and we envision that an analysis of the extended broadband emission spectrum such 
as that described in this paper will be similarly fruitful in terms of constraining the properties 
of the planetary atmosphere. Further EPOXI targets for which we will provide constraints on 
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the geometric albedo are WASP-3, TrES-3 and HAT-P-4, all of which will have Spitzer IRAC 
measurements of the infrared secondary eclipse depths in at least the 3.6 and 4.5 micron channels. 
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